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Preface 


We present here some of our recent work on structural identification and large- 
angle maneuvers with vibration suppression. Our recent work has sought to balance 
structural and controls analysis activities by involving the analysts directly in the 
validation and experimental aspects of the research. We have successfully imple- 
mented some new sensing^ actuation^ system identification^ and control concepts. An 
overview of these results is given herein# 
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we distinguish between modal and model identification. In modal identification, 
we seek to extract the natural frequencies, damping factors, and mode shapes from 
(typically) a free decay response to impulsive initial conditions. In model identic 
fication, we seek the estimates of a particular set of physical parameters contained 
in our best available model of the structure, so that our computed (free or forced) 
response of the system is in best agreement with the measured response. 

In our approach to model identification (see next page for roadmap figure), we 
utilize modal identification as a preprocessor. That is, we consider the 
{(io’s, and (j>*s) determined from free vibration measurements as additional 

measurements to our model identification process wherein we seek to fit the 
measured input/output behavior by the frequency response from our model of the 
systemc See reference 1 for details. 
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A Novel Approach to Structural Identification 


STRUCTURAL MODFX IDENTIFICATION 

Esiimaiion of the bcsi-fitting linear structural 
model via Creamcr/Junkins substruciuring 



o 


I FREQUENCY RESPONSE IDENTIFICATION 

Estimation of the frequency response function 
via Creamer’s algorithm 






MODAL IDENTIFICATION 

Determination of natural frequencies & mode 
shapes via Juang/Pappa Eigcnstruciurc Reali- 
zation Algorithm (ERA) 



INERTIAL TARGET TRAJECTORIES 

Compute inertial time history of each target’s 
coordinates via stereo ray intersection 




FOCAL PLANE TRAJECTORIES 

Connect successive image coordinates to form 
focal plane trajectories, for each of two (or 
more) focal planes 



best Jilting linear system: 

M X + C X + Kx = Bu 



measured natural 
frequencies & mode shapes 

CO,- ()),■ , i = 1,2, ... , 11 


field of view 1 field of view 2 











“ V. 




o 

IMAGE CENTROID CALCULATION 

From the synchronized cameras’ digitized 
image boundaries, determine die centroids of 
each image 

VIDEO PROCESSING 

Locate and digitize the image boundaries via 
Motion Analysis' edge detection processor 



camera 1 









This is a sketch of a 5 ft by 5 ft aluminum grid we*ve been using for a struc- 
tural identification and vibration control experimental test article* The grid is 
cantilevered in the vertical plane* We cut this grid from a single sheet of alumi- 
num, so it is jointless* However, a lap joint grid of identical geometry has been 
developed and is the subject of a similar experimental program by Alok Das et al. at 
AFAL. 

Note that we have 20 optical targets on the grid. Presently these targets are 
actively illuminated by fiber optic light guides, but we have also successfully used 
passive targets made from reflective tape ( "scotchli te” made by 3 m). These optical 
targets are visible in two video cameras (200 or 60 frames per second); stereo tri- 
angulation of centroided image coordinates leads to measured inertial trajectories 
for each target. We also have six strain gauges mounted on the back of the grid at 
the stations shown, and we utilize three piezo-electric accelerometers that can be 
mounted at any three of the grid locations* 

Excitation is provided by three grid-mounted reaction wheels (driven by Clifton 
Precision motors) with their torque axis in the plane of the g^rid in the directions 
shown (±20 oz-in* over the bandwidth from 0 to 40 Hz)* We also utilize an impulse 
hammer to impart calibrated initial impulse excitation at any station. 




We summarize the geometric equations governing triangulation of video camera 
imagery. With reference to the left figure, the image coordinates (x,y), measured in 
the positive focal plane of a single lens camera, are related (ref. 2) to the object 
space coordinates (X,Y,Z) of the imaged point, the camera's principal point 
(XcrYcfZc)f orientation angles principal point offset (xQ,yQ), and focal 

length (f), by the colinearity equations 


r c„(x-x.) + c,,(y-n) + c„(z-z,)i 
+ c„(>-n) + c„(z-z,)J 

= .n, j,./) 

”Cn C|| Ci)"l 
Cji Cn Cn I = 

poo-. 

= 1 0 cosy siny I 

r~cos0 0 
0 10 

ncos^ sin^ 0“| 
1 cosi 0 1 

I” Cn (X~Xe } + Cn(K— Kc ) + Co (Z~Ze )~l 

L ) + C„(i'-K ) + C„(Z-ZU 

j = G(X,r^.X,.K,Z, O', ./) 

^Cji Cn CnJ 

Lo -siny cosyJ 

La/i9 0 cos9~i 

L 0 0 J 


(1 ) 


We adopt a double subscript notation for eqs. (1) to denote the image 
nates of the i^^ point measured in the j camera's image space as 


coordi- 


^ojryojyfj) ^ 

yif — CJ (X| f Yf f Zi*, Xfy , Vcj , , <|>y , 0y , t x^j ♦ yc^ t fj) 


In the event that the camera position, orientation, and calibration coristants 
are considered known, measurement of image coordinates of object space points at 
unknown locations, eqs. (2) provide four equations" for the three unknown coordinates 
of each point. Equations (2) can be inverted by least squares to determine the ob- 
ject space coordinates. Dynamic triangulation must be preceded by a static calibra- 
tion to determine the camera calibration constants, and for the case of fixed cameras, 
the position coordinates of the cameras. The calibration should use at least three 
fixed targets that will subsequently be visible in the dynamic experiments. These 
three points serve to define the object space coordinate system. The points are 
numbered in ^ome arbitrary fashion. We use Point 1 as an arbitrary origin: 

(X^,Y-|,Z-j) = (0,0,0). Point 2 is used to define the X axis: (X 2 ,Y 2 /Z 2 ) = (X2r0,0). 
Point 3 is used to define the X,Y plane: (X 3 ,Y 3 ,Z 3 ) = (X3,Y3,0). Thus the first 
three points have a total of three unknowns (X 2 /X 3 ,Y 3 ), and each additional point 
introduces three unknowns (Xj^ ,Yj^,Zj^) for a total of 3N - 6 unknown object space 
coordinates. Notice that each object space point has four associated measurements 
(two measured coordinates in each of the two image planes). Thus, in the most gen- 
eral case, we have the 3N - 6 unknown object space coordinates plus the 18 unknowns 
associated with the cameras , Y^^ ,Z^^ ; for 1 = 1,2), We 

conclude that we have a total of 4N eqns. and 3N + 1 2 unknowns; if N > 12 conjugate 
images are measured, we have enough equations to determine all of the 3N object space 
coordinates and the 18 camera position, orientation, and calibration parameters. 
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Also of significance, when using Charged Coupled Device (CCD) focal plane arrays (as 
in the present application), the effective focal length is different for the x and y 
axes (actually, there is a ratio of the x and y dimensions of the array as a conse- 
quence of the rectangular pixels, which can be absorbed into effective focal lengths 
to be determined during calibration) thereby increasing the number of unknowns by one 
for each camera and one additional object space point is required. We have found 
this calibration process to be very well behaved and can be routinely accomplished in 
15 minutes of real time. 


The calibration process is strengthened, of course, by making redundant measure- 
ments and using other a priori measurements of object space points and/or camera 
calibration parameters. The precision of the results is dependent most heavily upon 
the accuracy with which the image coordinates are measured and, of course, the geo- 
metric strength of the triangulation process. Upon completing a satisfactory static 
calibration, the subsequent dynamic triangulation process need consider only four of 
eqs. (2) at a time to sequentially accomplish the least-squares solution for the 
object space coordinates of the points imaged on the flexible structure. 

The major hardware elements are sketched below and described briefly in 
the table. 
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Here we show two typical measured strain histories at two stations in response 
to impulsive excitation via the impulse hammer. The time span of these records is 
2.5 seconds; the sample rate was 700 Hz. An initial impulse was apj^ied" at“ the“ lower 
left corner of the grid. Notice from the fast Fourier transform (FFT) of strain, __ 
gauge #4 that the response is dominated by the first two modes, but significant exci- 
tation was imparted to most of the modes below 40 Hz. On the other hand, strain 
gauge #1 has a large contribution by modes 3 and 4. This is not surprising, because 
it is intuitively clear that the first torsion mode (mode 2, wilJi nominally zero 
strain in the center member) even though excited has near zero strain at station #1, 
but the second cantilever mode (mode 3) and the first "bow” mode (mode 4) are heavily 
excited by an impulse at the lower left corner, and these modes have significant 
strain at station #1. Based upon these and similar experiments, we conclude that" the 
strain gauges provide excellent sensitivity to all modes below 50 Hz, for impulsive 
excitation, and simultaneously using strain measurements at the six stations shown 
provide excellent observability for all modes below 40 Hz. 

As evident in the next figure, however, we have encountered some signal- to- 
noise difficulties when using the strain gauges to measure response to the 
reaction wheel excltatlono ^ 


Response of the Grid to Impulsive Excitation 



[calibrated initial impulse applied at lower left corner of the grid] 




Below we show two typical measured strain histories at two stations in response 
to harmonic excitation by reaction wheel #1, A 2-Hz harmonic torque is applied to 
the structure at the center of the bottom member, about a nominally vertical axis. 

The torque was applied for approximately 2 minutes to allow a steady state to be 
achieved before we acquired the above 2.5 seconds of data. 

Our experiments indicated that the Clifton Precision motors deliver very clean 
harmonic excitation, so the high-frequency variations evident in the strain-gauge 
response are virtually all due to noise. Strain gauge #4 gives data with a satisfac- 
tory signal-to-noise ratio. It is evident however that the much smaller strain 
levels and their closer proximity to the electrical disturbance of the actuator re- 
sulted in much noisier data from strain gauge #1. In both cases, the physically 
dominant steady-state response at 2.0 Hz obviously dominates the strain measurements. 
However, it is obvious that the "noise modes" are almost negligible in strain gauge 
#4, whereas they are very significant in strain gauge #1 output. Note the cluster at 
60 Hz in the FFT of strain gauge #1’s response. This is a near-certain indication 
that electrical noise from the motor is generating a significant part of the problem. 
We anticipate that better grounding of the motor will reduce this noise source, but 
there is the difficulty that locally small strains (which are associated with certain 
motions) are the fundamental cause of the poor signal-to-noise ratio when using 
strain gauges as distributed vibration sensors. 

These data indicate that the strain gauge is potentially useful in the present 
application, but other sensing approaches should be explored to eliminate signal-to- 
noise difficulties. As is evident from the impulsive response in the figure below, 
the optical sensing methods yield very clean deflection measurements of the low- 
frequency vibration (<20 Hz), and are therefore very well suited to the present ap- 
plication. We also show below a typical accelerometer record and its FFT. We have 
found the accelerometer data accurate between 2 Hz and 100 Hz, but very poor below 
1 Hz. Our present optical system and triangulation deflection measurements are not 
compatible with real time. Thus the optical system as we have it presently con- 
figured and implemented is suitable for system identification but not feedback 
control measurements. 
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Volts M Displacement 


Comparison of Video-Derived Position 
Measurements with Accelerometer Measurements 


1.5 


Video Data 
J9-1 



0 


Time ( sec ) 


10 


FFT of Video Data 
J9-1 



Frequency (Hz) 


Accelerometer Data 



Time (sec) 


FFT of Accelerometer Data 
Test J9-1 



Frequency (Hz) 


The graphs shown above provide some insight into the advantages and disadvantages of the camera 
system versus an accelerometer. The video data is very useful at the low frequency large amplitude end of 
the spectrum. The accelerometer response is most reliable at the high frequency end of the spectrum. 
This data was taken on the AFAL structural identification test article which is similar to the TAMU 
structure. The video data was taken with a TAMU RCA TC2811 60 Hz video camera. An AFAL En- 
devco model 7751-500 accelerometer was also used. 
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Our structural identification approach is based upon the following observations: 

• The input/output behavior is uniquely captured by the system transfer function 
(or frequency response function) matrix^ whereas an infinity of structural 
models has the same transfer function. 

• The most important necessary condition of a good model (for control purposes) is 
that we accurately model the actual system’s frequency response over the fre- 
quency range of interest. 

• Of all the models (realizations) we might use, the most comfortable approach is 
to modify in some "intelligent” way the model that grew out of our modeling 
effort (e.g,, a finite-element model of the structure). 

It is desirable that both free response and forced response behavior be 
accurately captured by the identified model. 


Basic System Realization Concepts 

Linear Aulononious System 

Notes 

X = Ax + Bii 


y = Cx 

The triple (A,B,C) is not unique. 

Time Response 

But G{s) is unique. 

x(/) = x(/„)t + 


y{t) = Ce'"'' x{to)[ + llCe^^‘'^Uhi{x)dx 

(A, ,B,,Ci) and {A 2 , /? 2 , C 2 ) arc said to be 

Laplace Transform 

equivalent if any of the following statements 
arc true; 

ji'(i) = {si - A y’ Bu{s) 

• The transfer functions arc equal; 

y{s) = Cx{s) = C{sI-Ay‘Bu{s) 

G i{s) = G2{s), lor aUs 

y{s) = G{s)u{s) 

• The weighting patterns arc the same: 

Transfer Function 

C//'7i/ = C2c'^^'B2 
• CiA\Bj = C 2 A* /^ 2 ) for all A:. 

G{s) = C{sI -A)' B 


Frcq. Response Function 

=> Any {A,B,C} which produces the correct 
transfer function C(5) - C{sl - A)‘ B is 

G (/(o) = C(jo)l - A ) ' B 

said to be a "realization" of the system. 
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Consider the case that a force or moinent is applied at a pointy but the response 
at many measurement stations is available* For this case, B and therefore G(ju) are 
vectors* Suppose that the frequency response function (vector) is measured over a 
frequency range ^ ^ ^ ^ax^ frequencies r 

notation: 


G = 



G((Oi,p)^ 


> = measured FRF, G(p) 




L = computed FRF using model vector p 


(3) 


and 


Q = {Oi ^2 ••• Y ~ measured free vibration Q’s, ^2 •••} ~ computed Q s using/;. 

(4) 


we seek the optimal estimate of the model parameter vector p which minimizes 


} = i j AG^ {p,(a)WGi<o)AG {p,(o) d(o + , AG(p,(o)=G((o) - G(o),p), Afi= Q - QQf) 

(5) 


For FRF measurements available at discrete frequencies, the integral can be re- 
placed by a discrete summation; we seek to minimize a weighted sum square of the 
residuals between all measured and modeled FRF's and fl's: 


AG = 


G{ta,)-G{(S>,,p) 
|G(oiu)-G((Ou,p) : 


FRF residuals, 


An = 


n, -n,(r)' 

' tij -Oitr)' 


free vib. frequency residuals 


( 6 ) 


and we are led to the least-squares differential correction algorithm: 




+ ... 


Ap 


Iaq 


-dG- 

[. where A= , 

J L^J 


and Pnew ~ /^old AP 


(7) 


Potential trouble. This approach works great, but only if: (i) the model is "good" 

and (ii) p ^ .is "close" to p. 
start 


375 



Creamer’s Tlii ee Step Identification Process 


Step 1: Identify a subset of the measured frequencies and motle shapes wiiich correlate 
well with the corresponding modeled frequencies and mode sliapcs, do a least square 
correction correction of the stiffness paramclcrizalion to improve correlation if 
necessary. 


Slep 2: Find the normalization of the cigcnvcclors (mode shapes) which results in the 
best least square fit for the FRF, i. c., find [gq, Gj, O 2 , . . . <'?,„+/} lo minimize 

.1 V^P‘l 0 “ ^^^PPmoilcI 0 j J 01 S ijfpcj (/ CO^ ) — (/ CO^ )j 


where 


^^P‘lmoM (/®) 




73 + ^ Or +a„t+\ 

CO '=1 ~ ..2 

CO,. -CO 


Slep 3: Estimate values for the linear mass and stiffness parameterization lo satisfy the 
orthonormality conditions in a least square sense; this leads lo the following pair of 
linear equations for [ir,Gii(l tc, in M ~ + S,- p,M,- , K = K^, + Z, ic,/C, : 


Here we show the first five free vibration mode shapes. In the table, we list 
the a priori finite-element model prediction of the first five natural frequencies, 
the values recovered from using the Ei gens true ture Realization Algorithm to reduce an 
impulse response of the grid as measured by the stereo triangulation process. As is 
evident, the agreement is generally good with the most significant prediction errors 
in the higher modes as might be expected. Applying substructure scale factors to the 
members contribution to the mass and stiffness matrices resulted in an identified 
structure with near-identical measured and modeled natural frequencies and mode 
shapes. The measured modes differed most significantly in that the first mode was 
observed to have a slight degree of torsion whereas the modeled first mode was a pure 
cantilever mode. The discrepancy was traced to a slight bend near the upper right of 
the structure, apparently this occurred when the structure was being mounted. 

Using only the optical data, we have demonstrated that the free response and 
input/output behavior of this structure can be accurately identified using the stereo 
triangulation system and the methodology we have developed. The presently imple- 
mented system works well only for frequencies below approximately 15 Hz. The dif- 
ficulty at higher frequencies is not a consequence of the camera frame speed 
(200 frames/sec), but rather because the limited spatial resolution of our cameras 
cannot "see” the small amplitude vibrations associated with the higher modes. 
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The appendage/hub structure was designed, fabricated, and assembled so as to 
produce a structure that is symmetric with respect to physical and geometric param- 
eters and that is also subject to a minimal dynamic asymmetry. Phasing problems that 
exist between appendages are apparent only after free vibrations are allowed to con- 
tinue for several minutes. This permits free and forced vibration data to be re- 
corded for a structure that exhibits an initial dynamic symmetry. 


TAMU Maneuver and Vibration Suppression Experiment 


PRIMARY SENSORS 

TORQUE SENSOR 
Sensor Developments 
custom unit with 
siffnal conditioner 


STRAIN GAGES 

Micromeasurements Div 
Measurements Group 
Std. metal film 


ANGLE ENCODER 
Teledvne Gurley 8700 
3G000 counts / rev 
C-TEK counter 
mod. LIN-101-5021410E 



^ECONDaS' SENSORS 

TACHOMETER 
ON REACTION WHEEL 
PMI Motion Technologies 


\ 

MOTOR CURRENT MONITOR 

hard wired to 
A/D board 


MOTOR VOLTAGE MONITOR 

hard wired to 
A/D board 




The below maneuver is from a student (Robert Byers*) recent M, S* thesis 
(ref, 3), Byers implemented a variable-structure feedback control law that roughly 
approximates a smoothed bang-bang near-minimum-time maneuver, but also exhibits good 
vibration suppression characteristics as is evident* The variable-structure approach 
is one of several robust control concepts we have studied in recent months, and we 
are engaged in an experimental effort to parallel and support our analytical 
activities* 


As an alternative approach, we show below an output feedback control law that is 
based upon a PDE description of the system dynamics, and we prove global stability 
using Lyapunov’s second method* This method is found to be very robust and, since it 
does not require a state estimator and since no spatial discretization approximations 
are introduced, it does not suffer from spillover; we feel this is a very attractive 
approach. We have shown that a near-minimum-time torque-shaped open-loop maneuver 
can be introduced and the Liapunov based method leads to a tracking-type feedback 
control that is theoretically globally stable and is, in fact, very robust with 
respect to modeling and implementation errors* It is this latter idea upon which we 
have based our early experiments; we will therefore go over the key ideas underlying 
this approach. 


Icxas A&M Maneuverable Space Structure Experiment 


15° Resl-fco-Rest: Maneuveu' 

Ti = .35 sec, a - .93, ft = .97, ij) ^ ip ^ 
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The below figure shows the experimentally achieved hub rotation angle as a 
function of time for a 40° near-minimum time maneuver. The first half of the maneu-_ 
ver was almost identical to the corresponding simulation (slight lag developed due to 
a lag in the compensator circuit of our power supply). However a large transient 
occurred at precisely the instant (maneuver mid— point) the commanded current reversed 
to initiate the breaking phase of the maneuver. This is because our low-budget power 
supply was unable to accurately generate the commanded current in the face of the 
variable load presented by the motor near torque reversal. During the last half of 
the maneuver^ the power supply compensation circuits and control commands combined to 
arrest rotation and vibration with a moderate overshoot. 

While the problems with our power supply are evident, we feel pleased with this 
first implementation of our approach and anticipate substantially improved results in 
future experiments. 


A Large Angle Maneuver of the TAMU Flex 

(Test Run No. 3 of 01/25/89) 
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Concluding Remarks 


• Novel Vibration Sensing and System Identification 

Methods have been Developed and Demonstrated 

• Large Angle, Near-Minimum-Time, Feedback Control 

for Flexible Body Maneuvers have been Developed and 
Successfully Demonstrated 

• Both of the Above are in a Preliminary State of 

Development, but Our Results to Date are Significant 
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